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are also looked at with the aim to determine their reliability and, eventually, 
why they don't work. Further aspects relative to the failure of the Operator 
Expansion Method to reproduce the results of more elaborate calculations are 
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I. INTRODUCTION 



The importance of the double beta decay process is well recognized. First, the neutri- 
noless mode, yet unobserved, is of fundamental interest, as it will be a signal for neutrino 
mass and lepton number non- conservation. Second, the double beta decay with two-neutrino 
emission ((3(321/), allowed in the standard model, is a very rare process which has not been 
experimentally observed until 1987, in the decay of 82 Se studied by Elliot et al. ffl. Subse- 
quent results in other nuclei were obtained by other groups (see || for a recent review of 
the experimental situation). Recently Balysh et al. M have measured the double beta decay 
half-life of 48 Ca. This nucleus is the lightest one for which such a measurement is feasible. 

On the theoretical side (see {||| for a recent review) the /3/32z/-decay which, at the 
beginning, was a well defined process in the standard model, has revealed as a real challenge 
for nuclear model practitioners. There are two reasons to this situation. On the one hand, 
the decay mode is highly suppressed and sensitively depends on poorly determined parts 
of the nuclear interaction. On the other, it is a second order process, which implies a 
summation on intermediate, and not always well determined, states. Thus, even if it is not 
a process involving new fundamental physics, the (3(32v-dec&y is related to a new type of 
nuclear matrix element. This one incorporates information on the wave functions that is 
not given by other standard observables. 

The difficulties in the calculation of the (3 /32z/-decay have been expressed in several dif- 
ferent but related ways in the literature: 

• In QRPA calculations, it has been related to the extreme sensitivity of the transition 
matrix element to the so-called g pp parameter, which governs the pn-excitations |(| § . 

• In the SU(4) language, it has been connected to poor determination of the nuclear 
force in the L = 0, S=1,T = 0(J=1 + T = 0) channel |9|f|. At the same time, it 
was also observed in this scheme that strong truncation of the model basis can produce 
undesirable contributions to the transition matrix element. 
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• It has been also connected to the bad description of the f3 + decays, which processes 
have been used to fit the unknown parts of the nuclear force [[7]. 

In order to avoid some of the previous uncertainties, an alternative approach has been 



proposed, the Operator Expansion Method (OEM) |T0|-|I"2|. 

In the present paper, we will focus our attention on the role of different parts of the 
nuclear force in the (5f32v-decay as well as on different methods used in the literature to 
describe this process. For that purpose, the simplest nuclear transition to be studied is the 
/3/?2z/-decay in 48 Ca, that offers a double advantage. There exist both a sensitive experi- 
mental value (Tfy = U.3tH [Stat] ± 1.4 [Syst]) x 10 19 y [|) and an elaborate shell model 



calculation [13], which is the natural calculation scheme for this nucleus. While doing these 
studies, we will have in mind a long standing problem. Different calculations were approx- 
imately leading to the same decay rate whereas the intrinsic sign of the transition matrix 



element was not the same jnj, requiring some clarification. With this respect, we will in 
particular show how higher order effects in the SU (4) symmetry breaking interaction modify 
previous estimates. On the other hand, a critical study of the OEM approach has been done 
in ||15|| . Based on an analysis of the SU(4:) symmetry breaking effects, other features of the 
OEM have been revealed, which deserve discussion. Our work is therefore concerned more 
with the role of various approximations than with a realistic calculation of the process. We 
will use an analytical force to achieve this objective. This allows us to easily compare various 
methods of calculation and to switch on and off the different parts of the interaction. 

The plan of the paper is as follows. We remind in the second section expressions for the 
standard transition operator and that one obtained in the OEM approach. In section 3, we 
introduce in the OEM expression the Coulomb splitting effect and make a comparison of the 
corresponding result with the transition operator derived independently at the first order 
in the SU(A) symmetry breaking interaction. The effective NN interaction that we will use 
for our study is specified in section 4. The fifth section is devoted to a presentation of our 
results together with a discussion. 
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II. THE /3/32i/-DECAY: DESCRIPTION OF THE TRANSITION OPERATOR 



The /3/32z/-decay is allowed in the Standard Model. In this process, two neutrons decay 
to two protons with the emission of two electrons and two neutrinos. The Lagrangian 
responsible for that process is the standard Fermi one. Under the usual assumptions (the 
impulse approximation is assumed, the lepton energies are replaced by their average values 
and non- conservation of isospin is discarded), we obtain that the life time is expressed by 



Tf/fJ ' = G 2v \M gT \\ (1) 

where G21,, which can be found in contains all the leptonic part and the integral on 

the phase space. The nuclear information is included in the nuclear matrix element 
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Mar = - I dte lAt < 
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iHt Ae- iH \A 



10+ >, (2) 



with A = J2&i T ti ^e Gamow- Teller operator, and A = ~ (Ei — Ef). 

Methods for evaluating the /?/32^-decays differ in the approximations made in order to 
calculate (||]). What we will call the standard method consists to insert a complete set of 
intermediate states in the commutator This allows one to perform the time integral and 
we obtain 

„. _ < OtU|l+ >< 1+|I|0+ > 

Theoretically this method is exact, but in practical calculations some limitation occurs. 
Estimating the matrix element given by Eq. (|3|) implies to consider all the intermediate 
states contributing to the sum. Thus, the method could not be the most interesting one 
if the intermediate states or their energies cannot be well determined and if cancellations 
between different contributions are present. As we will explain later, the approximated 
SU (4) symmetry of the nuclear forces tells us that such cancellations must be present in 
(0), leading to uncertainties in employing this method of calculation. 

On the other hand, in general, nuclei which undergo double beta decay are open shell 
nuclei and the usual formalism for describing them is the QRPA. In that case, there are 



new difficulties in the evaluation of the Mg\ amplitude due to the fact that the QRPA is 
near the collapse for the physical values of the so-called g pp parameter. Improvements of the 
QRPA, like full-QRPA, renormalized QRPA or full-RQRPA, can solve these problems but 
other ones arrive, like the violation of the Ikeda sum rule (see discussion in ref. @||, and 
see also flT7|] ). 

As mentioned in the introduction, the above difficulties have been a source of concern. 
An alternative approach to calculate directly the exponentials appearing in @, the operator 
expansion method (OEM) |]TQ| |l2j| , has thus been proposed some years ago. This is impossible 



in a general manner and the main assumption made by the authors is to only retain two- 
body operators, which supposes that only the two nucleons involved in the transition are of 
special interest in the calculation. This is like a spectator approximation in the sense that all 
the nucleons not involved in the (5(3 transition don't contribute to the process. On the other 
hand, the interaction between the two active nucleons is included to all orders, neglecting 
some parts as we will explain later. A diagrammatic view of this approach is shown in Fig. 
0a. 

In the simplest approximation, the kinetic energy term, the spin-orbit and the tensor 
potentials are neglected in the Hamiltonian H appearing in (^j. Under these assumptions, 
only the central part of the potential contributes. Starting from its expression written as 
follows 

H ~ V c = [v (r) + v T (r) n.fj + v a (r) a^Bj + v aT (r) 7*. 73 a^Bj] , (4) 

i>j 



the /?/?2z/-transition amplitude in the OEM approach, first derived by Simkovic et al. [1C 
and Ching et al. fill , can be expressed as 

M e °# M =<0+|EM li r+r/|0+>, (5) 



i>3 



with 



24 [v a (r) - v T (r)] 

Ma = — — ill) 

3 A 2 — 16 [v a (r) — v T (r)] QKJ} 

+ § [2rw (r) - y a (r) - v T (r)] 

A 2 - 16 \2v aT (r) - v a (r) - v T (r)] 2 KJJ ' K ' 



where Qo,i represents the projector operator on spin and 1 subspaces. 

The Gamow- Teller operator is a SU(4) generator. Then, from equation (Q), we observe 
that Mgr will be zero if the Hamiltonian were SU(4) invariant (we discard transitions 
between members of the same SU (4) multiplets). This result is not obvious from (|3|) where, 
in the general case, each intermediate state could give a non-zero contribution, the zero being 
obtained after summation over all intermediate states. Usual nuclear potentials are not so 
far from the SU (4) symmetry, then we must expect cancellations in the summation present 
in ([H). From this point of view, expression (g|) is more transparent, because the condition 
over the central force to be SU(4) symmetric is just v a (r) = v T (r) = v aT (r). Notice that 
in a particular case, actually close to most realistic transitions, the vanishing of the matrix 
element would result from the fact that the operator A in Eq. ([3D acting on the final state 
gives zero. 

III. OEM AND SU(4) APPROXIMATION 

We can improve the OEM model, Eq. (|6|), by incorporating the contribution of the 
Coulomb interaction and then make a comparison of the expression so obtained with that 
one derived in the first order SU(4) symmetry breaking approximation. 

The Coulomb interaction is an important ingredient of the spectroscopy of medium- 
heavy nuclei involving different charges as it provides a few MeV shifts, which compare to 
the energy splittings produced by the strong interaction itself. For our purpose, we added 
to this interaction, H s , a constant term proportional to the third component of the total 
isospin of the nucleus 

H = H S + A C T 3 . (7) 
The above Coulomb force can be easily included in the operator appearing in (||). Using 
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that 



T 3 , A = A, we obtain the modification of (|6]) into 

24 [v a (r) - v T (r)] 



7vf 



J (A + A c ) 2 -16K(r)- Vr (r)] 2 

8 [2?y r (r) - (r) - v T (r)] Q 
' (A + A c ) 2 - 16 [2 W(7T (r) - u CT (r) - v T (r)] 2 1 



(8) 



From the definition of A and @, we have the relation, A + A c = | [E* — E S A , where Ef is 
the strong interaction contribution to the energy of the state. 

As we said before, the nuclear forces are not far from the SU(4) symmetry and, if that 
symmetry were exact, the double beta transition amplitude will be zero when it connects 
states belonging to different £{7(4) multiplets. In order to look at the consistency of our 
results, it can be useful to study the first order correction of our expressions in the 577(4) 
breaking parts of the force. Let us write for that 



H = H + H 1 , 



(9) 



where H (Hi) represents the 5^7 (4) symmetric (breaking) part of the force. The hamiltonian 
Hq is a purely central force and has two terms in the spin-isospin space, one proportional to 
1 and the other proportional to the Casimir of the £77(4) group, (fl.Tj + <7j.<7 j + fi.Tj (Ji.&j) 



(10) 



The Hamiltonian, Hi, which is able to give a double beta transition at the first order in 
the SU (4) symmetry, involves two other combinations of the components appearing in the 
central force, Eq. (01), 



+ 
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To 



( 2 

(v T (r) + v a (r) - 2v aT (r)) f Ti.Tj + Ui.&j - -Ti.Tj Oi.o, 



(11) 



Starting with (fj), we observe that there are two different situations. First, let us consider a 
nucleus with only two active nucleons. In that case it is obvious that [Hq, Hi] = 0, and the 
exponential of the Hamiltonian present in (H) can be splitted in two exponentials relative to 



Hq and H\ respectively. The first one, related to H$, commutes with all the operators and 
can be ruled out. The second one, related to Hi, can be expanded to the first order in H 1 
and we obtain 



Mi 



(2neutrons) 



gr 



f <o+| 



2 (A + A 

Calculating explicitly these commutators, we get 

1 



Hi, A 



A 



10+ > . 



M, 



(Ineutrons) 



gr 



(A + A c )' 



<0+|E{24 K (r)-v T (r)]fi„ (ij) 



i>j 



(12) 



+ 8 [2v aT (r) - v a (r) - v T (r)} Qi (ij)} t+t+ |0+ > . 



(13) 



This result agrees with the Coulomb corrected OEM expression, Eq. (8), when this one is 
expanded up to the first order in the SU (4) symmetry breaking. 

Surprisingly, a different result is obtained when we consider a nucleus with more than 
two valence nucleons. In view of its importance for a comparison with the above OEM 
result, Eq. (8), we give here some detail on its derivation. 

Beyond the two valence nucleon case, [H Q , Hi] doesn't vanish and from (||) it can be 
shown that, up to first order in Hi, 



gr 



J o ^ e ^+A c )^liJ_ <0 +| [£^[^[[ jF f 1) l] ; A]]...]|0 + >. (14) 

™ _1 n-l 



In order to perform the sum present in Eq. (0), let us define an operator E = Hi,A,A 
and introduce an operator C solution of the equation B = [H , C] . In terms of the C operator, 
(|i"4l) can be rewritten as 

00 [it) n 



dte i{A+A c )t <0 +| e iH t Ce -iH t_ C 

2 Jo f 1 L 

1 < + |C|0+ > . 



|o + > 



A + A c 



(15) 
(16) 
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The solution of the equation B = [H ,C] is 

POO 

C = -i\wy dq e iHoq B e~ iBm e~ Eq . 

£^0 Jo 

Performing this last integration and using the explicit expression of B, we obtain: 



(17) 



A 



10+ > . 



(18) 



2 (A + A c 

This equation can also be obtained from (^). To do that, we must realize that up to first 
order in the Hi Hamiltonian, the double beta decay implies a transition between different 
SU (4) multiplet states. In the case of interest here, it involves the state | [4, 4] T — 4, S — 0) 
associated with the 48 Ca ground state and the state | [2, 2] T = 2, S = 0) related to the 48 Ti 
0. Then, there is only one intermediate state contributing to the amplitude, which is the 
Gamow- Teller resonance of the initial state (with an energy E n = E t + A c ) : 

< 0+UI1+ >< 1+UI0+ > 1 



E 



< 0J\A.A\0i >, 



(19) 



„ {Et + A c ) - | (Et + E f ) ~ (A + A c ) " f 
and Hi appears in the mixing in the final nucleus between the two SU(i) representations 



10/ >= 10/ >o+E 



1 



|0+ > < o+i^!|o;> 



0- 



(20) 



(E f - E r 

where |0 + > is the pure SU (4) final state. When we introduce (p0|) in (|19D , only states 
|0 + >o belonging to the same SU(4) supermultiplet as \0f > (= \0f >o) can give non-zero 
contribution and the energies of these states are E r = Ei + 2A C . Then we obtain 

1 ^ 1 



E 



(A + A c ) V (E f - E r 

-5 < 0+|#iAA|0 + > . 



o< 0j\Hi\0 r t > < + |A.A|0 + > 



(21) 



2 (A + A c 

This result is in agreement with Eq. fllif) because the other terms present in the double 



commutator, < 



/ 



Hi, A ,A \0f >, vanish in the SU(4) limit. The main point here 



is that this expression has a sign opposite to (I2"|). This is due to the presence of many 
nucleons operators in the former expression, as we represent in Fig. |T]b, while the latter one 
only contains two-body operators with the consequence to provide the wrong sign in the 



first order SU(4) symmetry breaking limit. In particular, contributions due to pure Pauli 
antisymmetrization, as those depicted in Fig. |l]c, are not accounted for in the OEM. 

More important, in the simplest case where the operator Ho in Eq. can be approxi- 
mated by the sum of the single particle energies, the different commutators appearing in this 
expression can be calculated. Their contributions, which form a non- convergent geometrical 
series, are given, up to a factor, by the sum 



i(l + 2 + 4 + 8 + ...)(...), (22) 

where the first term in the parentheses is that one retained by the OEM. To get these 
contributions, we used the relation, E? — Ej = 2 ( A + A c ) . Formally, the above sum can be 
performed with the result 

—(■■■) = -(...)• (23) 
This is the result obtained from a direct calculation, Eq. flUf). It specifies in two ways the 



failure of the OEM demonstrated on a quantitative basis by Engel et al. [Hj) . i) Among the 



contributions that are accounted for by the expression, Eq. flT^D, it indicates which one is 
retained by the OEM. ii) The energy difference, Ef — Ej, implies the single particle energies 
of nucleons in the initial and final states. These ones involving the core particles, it makes 



it clear that the various commutators appearing in Eq. (jTj) involve three and more body 
operators. 



IV. THE EFFECTIVE NN POTENTIAL 

As mentioned in the introduction, we are motivated in this paper by two different points. 
First, we want to study the contributions of the different pieces of the nuclear force to the 
/3/32z/-decay. Second, we want to compare the two calculation methods presented in the 
previous section. We will focus on the double beta decay of 48 Ca because it is a nucleus 
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which can be theoretically described in the nuclear shell model and we can do reliable 
calculations with both methods. To accomplish our objective, we use an analytical force 
which, therefore, could not be the best one but, as we will observe later on, the results are 
good enough to make it credible. In this way, we can easily connect and disconnect the 
different pieces of the force and calculate the matrix elements of the operators present in 
(g). We have performed our calculations using the OXBASH code ( flSH). 

The shell model space is the full fp shell with the single particle energies e/ 7/2 = 0, 
e P3/2 = 2.1MeV, e Pl/2 = 3.9MeV and e/ 5/2 = 6.5MeV. We have used for our calculations 



the Bertsch-Hamamoto force |]19[ . This force has a central part which in momentum space 
is given by: 

f f 



2 r-, 2 
1 . . mi 



Vc (q) = \ — I 77<?1.<7 2 f x .f 2 " 2 + 

m n J [6 q z + m^ 

(a a n 5 =°n r=1 + ^n^n^ ) + 

v ' q l + m l a 



(a fe n 5=0 n T=1 + £> b n 5=1 n 



T=0 



m 



2 



q z + m^ 

with / = 0.97, m n the pion mass, m a = 2.5fm _1 , mj = 4fm _1 , and a tensor part 
Vt (q) = - [ — ] \ (3ai.qa 2 .q- o x .a 2 q 2 



(24) 



x 



/ 3 

T\.T 2 



tl \ to 



(25) 



^q 2 + m\ q 2 + m\ ) q 2 + m 2 , 

Bertsch and Hamamoto (B.H.) fitted the parameters of the force in order to reproduce Reid 
soft-core G matrix elements and used this force to estimate the Gamow- Teller strength at 
high excitation in 90 Zr. The values for these parameters are given in Table |. The main 
features of this force are: (i) it contains the one-pion exchange, which governs the long range 
part of the force, both in the central and the tensor terms; (ii) the other terms of the central 
force are pure S wave interaction; (iii) the attraction in the S — 1, T — channel is bigger 
than the one in the S — 0, T — 1 channel (thus the pairing in the T = channel will be 
greater than the usual pairing in the T = 1 channel); (iv) the B-H tensor force was fitted 
to be used only for L even waves (the T = channel) and in this order only one of to,i 
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parameters is necessary. The authors of ref. JT9[ chose to — while the tensor interaction 



in the T = 1 channel was completely discarded. We have checked that this tensor force, as 
an effective one, is consistent with the deuteron D-wave. In order to have a simultaneous 
description of the three nuclei involved in the transition, 48 Ca, 48 Sc and A8 Ti, we also fitted 
the average Coulomb displacement in Eq. (|7|), A c , using the relative position of the ground 
state of 48 Ca to that of 48 Ti. 

We nevertheless observe that the B-H force evidences some undesirable features when 
it is applied to the study of the nucleus spectroscopy. For instance, in the region of 48 Ca 
of interest here: (i) it gives a state density at low energy larger than obtained with other 

or 



standard potentials like modified versions of the Kuo-Brown G-matrix interaction pO 
PH ; (ii) the splitting between the first J = + , T = 1 state and the first J = 1 + , T = 
state of 42 Sc has the wrong sign; (iii) if we extend the tensor force as it is in the original B-H 
force to the T = 1 channel, it produces quite important matrix elements. What happens 
is that this force must be used with some short range correlations which will decrease its 
effective intensity. We will not introduce short range correlations and for this reason and 
from the fact that a so simple force cannot have unchangeable parameters in a large range of 
nuclei, we slightly modified the original B-H parameters, (i) we fitted a a and b a , reproducing 
the relative position of the ground states of 48 Ca, 48 Tz and the first J = 1 + state of 48 Sc. This 
change has reduced the state density in 48 Sc and has also corrected the splitting between the 
two first + states of 42 Ca as well as the splitting between the first + and the first 1 + states 
of 42 Sc (these states are important to determine the two-body effective interaction), (ii) we 
changed to and ti in such a way that the tensor matrix elements between two particles in 
the fp shell coupled to T = 1 has been strongly reduced but without change in the matrix 
elements of particles coupled to T = 0. These parameters are also given in Table | as modified 
B-H. Due to their limited number, our force is not the most realistic one. Results presented 
here cannot therefore compete with other ones which rely on a better force. As it can be 
seen from the results we obtain, they are realistic enough however so that our schematic 
study makes sense and can provide sensitive information. 
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In Table ||, we present the first states with quantum numbers 1 + for 48 Sc. Notice that, 
with respect to the state density argument presented above, this table is partly misleading. 
Higher energy states should be included in the comparison. 

We looked at the distribution of the Gamow- Teller strength for these forces and compared 
it with a standard calculation performed with a modified Kuo-Brown interaction [|20|| . As it 
can be observed in Fig. ||], there is no difference between the strength calculated with the 
B-H or the modified B-H interactions and that one using the potential of ref. [ 20| . The f3 + 



strength from the final state has also been looked at. It is shown in Fig. § for the same 
models as mentioned above. Its relevance has been mentioned several times in the literature 
and re-emphasized recently in ref. p^j . It represents an important constraint. We observe 



that for the B-H and modified B-H potentials the results are hardly distinguishable; for 
the modified K-B interaction more strength is concentrated in the low energy region. The 
essential point is that the [3 + strength is large where the Gamow- Teller strength is small and 
vice versa. Moreover, the contributions of the two regions could be opposite in sign, which is 
not observed in the strength, making difficult an accurate determination of the total matrix 
element. 

We want to stress here two points. Our modification of the Bertsch-Hamamoto force has 
not a fundamental origin. In that sense, we cannot say that this modified force is better than 
the original one, but it gives a better description of the spectra for the states of the nuclei we 
are considering. On the other hand, the main motivation for using the Bertsch-Hamamoto 
force, or a modified one, is that its analytical structure allows one to make a simple analysis 
of the role of the different pieces while the total transition matrix element is compatible 
with all previous calculations, as it will be shown later. 

V. NUCLEAR POTENTIAL AND /?/32i/-DECAY 

Our first study concerns the role of the central force in relation with the contribution 
of the spin-orbit. In this order, we turned off the tensor potential. The spin-orbit energies 
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are multiplied by a factor 7, running from to 1, in such a way that, for 7 = 0, only the 
central potential contributes while, for 7=1, the spin-orbit splitting energies are completely 
accounted for. The results for the different calculation methods (|3]), (0-8) and ( |i~8"|) are given 
in Fig. f§ 

Comparing Mg\ for the B-H and modified B-H potentials in absence of spin-orbit po- 
tential (7 = 0), we observe that small changes in the values of a a and b a make the double 
beta amplitude to go through zero. Hence, the contribution of the central potential by itself 
is not completely under control. This result points to the relative weight of the forces in 
the (S — 1, T — 0) and (S — 0, T — 1) channels, which plays an essential role in the present 
field. 

Quite generally, interaction models based on nucleon-nucleon scattering data have a 
strength in the (S = 1,T = 0) channel bigger than in the (S = 0,T = 1) one (as is the 



case of the B-H potential, see Table |HI[) , but most effective nuclear potentials fitted to 
reproduce the spectra of nuclei give a pairing for (S = 0, T = 1) states stronger than for 
(S = 1,T = 0) states. Typically, the situation is characterized by nuclear matrix elements 
like those displayed in Table JTJ, calculated for the B-H and the modified B-H potentials. 
It has a direct relationship to the relative weight of the forces, v T (r) and v a (r), in Eq. 
(U). The issue is an important one, which has a close relationship to the sensitivity to the 
so-called g pp parameter appearing in other approaches. As there, one has to hope that the 
fit of the effective nuclear potential model to a few relevant experimental informations will 
allow one to minimize uncertainties. In ref. p3j , Poves et al. considered the same problem 
in terms of two factors, A01 and A10 , multiplying respectively the strengths of the forces 
in the singlet and triplet spin channels. Starting from a force that was already good, the 
variation for these factors is actually smaller than what is suggested by the comparison of 
our matrix elements given in Table [TTT] for the B. H. and the modified B. H. forces. In ref. 



24[, one can find a recent study on the (pn) pairing and the relevance of the point here 



underlined in heavy nuclei, which are studied in the QRPA approach. This is also discussed 



in ref. [pq] . An argument is sometimes advocated for the change of the relative strength of 
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the forces in the singlet and triplet spin channels when going from infinite nuclear matter 
to finite nuclei. It relies on the effect of the spin-orbit force. The force in the singlet spin 
channel is coupling preferentially particles with the same quantum numbers, (j, I), whereas 
the force in the triplet spin channel rather couples spin-orbit partners, the effective force is 
favored by the absence of spin-orbit splitting in the first case while it is disfavored by its 
presence in the other. 

The role of the spin-orbit interaction has not received much attention in the field, proba- 
bly because it is known and is not considered as a free parameter. In an approach based on 
the SU (4) symmetry like that one referred to here, it has some relevance since it is a piece 
of the interaction that breaks the symmetry. Its importance can be seen by looking at the 
dependence of Mg\ on 7 shown in Fig. |j. It represents a quite important contribution for 
both potentials. In the B.H. case, it produces a change in sign while in the modified B. H. 
case, it enhances the amplitude by a factor 2. Algebraically, the effect is roughly the same 
and the difference by a factor 3 between the results for 7 = 1, Mg\ (B-H) = 0.044 MeV -1 
and MgJf (mod B-H) = 0.15 MeV -1 , is thus due to the central potential contribution, clearly 
over-estimated in the modified B-H potential. 

It is instructive to look at the detail of the contributions of the intermediate states to 
the matrix element Mgq-, in Eq. (|^). This is given in Fig. [| for the various models whose 
Gamow- Teller (3~ and (3 + strengths were shown in Figs. |2| and |3| respectively. For the B-H 
model (as well as the K-B model), there are contributions with both signs, respectively 
located at low and high energy. The dominant contribution in the low energy range is 
indirectly an effect of the spin-orbit interaction which brings down some states into this 
region and at the same time some strength. It is partly cancelled by a contribution in the 
Gamow- Teller resonance region, which is reminiscent of that one estimated in the <S77(4) 
symmetry approach (see below) or that one calculated in || on the basis of the dominance 
of this resonance in the sum entering Eq. (|3]). For the modified B-H model, all contributions 
are positive. The low energy range one has the same origin as above, whereas that one in 
the Gamow- Teller resonance region has the opposite sign. This is due to the change in the 
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relative strengths of the forces in the singlet and triplet spin channels evidenced by these 
models (Table [TTTD . As a result, the matrix element Mg%- for the modified B-H model is 
significantly larger. 

We now focus on the modified B-H potential and compare Mg% with Mj^ (4) . We observe 
that Mf r ^ gives a reasonable estimate of Mg^i up to a factor 2. But this nice result is 
partly due to the crossing of the two curves in Fig. |4], Mgj- and Mg^ 4 \ which makes their 
difference to remain in a relatively small range. The accidental character of the agreement 
is evidenced by looking at results for a different choice of the central force. Thus, for the 
case of the B-H potential, M|^ (4) varies from -0.072 MeV" 1 for 7 = to -0.053 MeV" 1 for 
7 = 1 while for the same potential Mg^- takes values from -0.014 MeV -1 for 7 = to 0.044 
MeV" 1 for 7 = 1. We must conclude that, even if mJ? (4) is a good estimate of the order of 
magnitude, it does not give the right sign and, moreover, differences for the absolute value 
can be as big as a factor 2 or more. The change in sign for Mg% , when going from the 
B-H to the modified B-H potential, is due to the relative value of the T — 1 and T = 
pairing as is shown in Table [TIT]. 

Looking at the dependence of Mg? on 7, we observe that its value is relatively stable. 
When 7 runs from to 1, the spin-orbit contribution to the wave function is included to all 
orders in Mg? but its vertex contribution, through the operator in @, is not considered. 
What shows the evolution of Mg^ is that this vertex contribution is the dominant one. 
Naively, we could conclude that a better estimate is to approach the Hamiltonian in (0) 
by H = H so + A C T 3 , but this contribution vanishes. In fact, the result we obtained for 
Mgtj- implies a large interference between the spin-orbit potential and the central potential. 
This is perhaps a consequence of the strengthening of the force in the singlet spin channel 
with respect to the triplet one, which we mentioned above as being indirectly due to the 
spin-orbit force. 

Looking; now at M°# M , we must conclude that its value is mostly independent of the 
spin-orbit potential. This statement is also true for the B-H potential. In that case, 
runs from -0.011 MeV" 1 for 7 = to -0.010 MeV" 1 for 7=1. The difference in sign for 
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Mg^ M calculated with the B-H and the modified B-H potential is again due to the relative 
value of the T — 1 and T = pairing. In this case, Mg^ M is not a good estimate; neither 
the sign nor the absolute value are well reproduced. In respect to the Mg^ M calculation, we 
must emphasize that the Coulomb force cannot be neglected. If we use (0) instead of (8) our 
results for M^ M vary from 0.0005 MeV" 1 for 7 = to 0.0003 MeV" 1 for 7=1, evidencing 
an absolutely non sense result. We observe that Mg® M and Mg^ 4 ^ have the same sign, in 
apparent contradiction with what was said in section 3. What happens is that Mg-f M has 
a peculiar behavior when the SU (4) breaking part of the central force is reduced, crossing 
the zero and changing sign when we multiply (11) by a factor k and study the limit of k 
going to zero. 

Beside the role of the spin-orbit potential, whose importance has been discussed above, 
we also considered the contribution of the tensor potential. We found that this one does 
not change the /3/32i/-transition amplitude in a significant way. Only a slight decrease was 
observed. This can be seen as due to an effective decrease of the spin-orbit interaction which 
is in fact observed around 48 Ca and has been attributed to the tensor force in the past [f2U . 
Our full results so obtained are: 

Mpj- = 0.135M6V" 1 ' (26) 

Mg¥ {4) = 0.084MeV-\ (27) 

M°f M = O.OllMeV" 1 , (28) 
for the modified B-H potential and 

Mf r = 0.032MeV-\ (29) 

M|^ (4) = -0.051MeV-\ (30) 

M^r M = -0.008MeV-\ (31) 
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for the B-H potential. 

Previous results for the amplitude calculated in the standard way are summarized in 
Table [TV| They can be compared to the experimental value Mgr = 0.074^°' 015 MeV -1 , also 
given in the Table. As it can be seen, our results are between a factor 2 too high for the 
modified B-H potential and a factor 2 too small for the original B-H potential. In view of 
the simplicity of the force used in present investigations, which has allowed us to study the 
role of its different pieces, results can be considered as reasonable. Summarizing the main 
features, it can be noticed that the amplitudes, M§ r and Mg^\ are shifted upwards by 
roughly the same amount when going from the B-H to the modified B-H potential. This 
is related to the change in the relative strengths of the force in the S = 0,T = 1 and 
S = 1,T = channels. The difference between these two amplitudes, indirectly due to the 
spin-orbit force, is relatively insensitive to this modification. Both effects are important to 
get a value that compares to the experimental one. The results for Mg^ M also show some 
sensitivity but are out of range in any case. Results of previous calculations employing other 
methods are given in Table V. It is seen that the above values fall in the range of the more 
realistic estimates, especially that one by Caurier et al. |T3(] , Mg\ = 0.065 MeV -1 , which is 
probably the most elaborate one. This indicates that our study, though schematic, deals 
with the real problems underlying the calculation of the Mgr amplitude. 



VI. CONCLUSIONS 

We considered the f3(32h / -decay process in the nucleus of 48 Ca with the aim to analyze 
different methods used in estimating the corresponding transition amplitude or to study 
the role of different components in the nuclear interaction. In this sense, 48 Ca is used as a 
theoretical laboratory for testing various approaches. In all cases, we have performed our 
calculations in the full f — p shell. In order to discuss the different terms of the potential, we 
used the B-H potential, which is analytical, and adapted it to our nuclei giving rise to what 
we called modified B-H potential. Our conclusions do not depend on particular aspects of 
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one or the other potential. 

First, we confirmed the strong sensitivity of the (3f32v amplitude to the relative strength 
of the potentials in channels (S — 1, T — 0) and (S — 0, T — 1) . This is well a known result 
and its importance has been emphasized in the QRPA calculations (in a different language, 
this was firstly pointed out by || as well as ||). But our main conclusion is that the 
central force alone, due to this cancellation, does not provide the leading contribution to be 
considered in the calculation of the (3)32 ^-transition amplitude. 

From our results, the single particle spin-orbit force appears to be the main ingredient 
in determining the actual value of the amplitude. Nevertheless, this term of the potential 
alone is not sufficient as the total amplitude then vanishes. Any sensitive estimate of the 
amplitude requires the interference between this spin-orbit term and the two-body parts of 
the strong potential. In this interference, operators involving three, or even more, nucleons 
could appear. 

Concerning the tensor potential, we did not find it was relevant for the j3 /32v-decay 
amplitude. 

The results obtained in the OEM approximation are far from the exact calculation. The 
OEM is not an approach under control, as already mentioned in the literature. Sizeable 
corrections come from three or more body operators but the non-convergent character of 
the expansion don't let much hope that the corrections are manageable. We have put in 
evidence that the OEM has a wrong SU (4) limit and this difference is also originated from 
the many body operators. Moreover, we observed that the Coulomb potential cannot be 
neglected at all in this scheme. 

For the future, one can imagine to improve the approach based on the SU (4) symmetry. 
It is not clear however whether accounting for the spin-orbit splitting is feasible while keeping 
a rather simple form for the expression of the transition amplitude. Another issue concerns 
the sign of the contributions of the different intermediate states to the total /3/32z/-transition 
amplitude. The (3~ and (3 + Gamow- Teller excitation from the initial and the final states 
respectively only know about the magnitude. The interesting question is to know whether 
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this information together with the knowledge about the total transition amplitude, Mgj-, 
can provide a clue as to the constructive or destructive character of the partial contributions, 
as exemplified by two estimates presented in this work. 

This work has been partially supported by DGESIC (Spain) under contract N°PB97- 
1401-C02-01. 
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TABLES 







a a 


h a 


CLb h t\ 


to 


A c (AfeL) 


Bertsch-Hamamoto 




-8.28 


-14.33 


6.56 11.20 0.89 







5.352 


Modified B-B 




-8.89 


-14.13 


6.56 11.20 1.335 


1.335 




5.512 


TABLE I. Parameters entering 


the nuclear force, equations (24, 25). 






Exp. 


6.51 


0.46 


0.54 0.64 0.74 


1.19 


1.48 


1.66 


B-B 


4.96 


0.27 


0.68 0.83 


1.17 


1.44 


1.67 


Modified B-B 


6.50 




0.55 0.93 


1.28 




1.60 


Modified K-B ref | 


20f 


6.49 




0.63 


0.97 


1.35 


1.67 


Modified K-B ref | 


21| 


6.38 




0.52 


1.07 


1.56 





TABLE II. In the first column, we give the position of the first state in 8 Sc relatively to the 
fundamental state in Ti. The following columns contain the excited energies of the following 1 + 
states in Sc relatively to the first 1 + state. 





S = 0T= 1 


S= 1T = 


B-B 


-2.98 


-4.87 


Modified B-B 


-5.03 


-4.24 



TABLE III. Diagonal matrix elements (f 2 L = 0, ST) for the cases (S = 0, T = 1) and 
(S = 1,T = 0) for the central part of the B-B and the modified B-B potentials, in MeV. The 
relative weight of this type of matrix elements determines the sign of the corresponding part of the 
force to the double beta amplitude 
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Mpj- (MeV- 1 ) 




M|f (MeV^ 1 ) 




Mpj- (MeV" 1 ) 


Mod. B-H 


0.135 


ref 


13| 


0.065 


ref [26| 


0.089 


B-H 


0.032 


ref 


15| 


0.043 


ref [27| 


0.15±0.07 


Experimental 


074 + - 012 

u - u ' -0.015 


ref 


20| 


0.055 







TABLE IV. Double beta amplitudes for the two v mode. The second row of the first column 
corresponds to the modified Bertsch-Hamamoto force, the third row of the first column to the 
original Bertsch-Hamamoto and the fourth row of the first column to the experimental result. All 
the other columns correspond to usual shell model calculations except for ref. |2?]] which corresponds 
to a shell model Monte Carlo technique. The experimental result corresponds to || . We have taken 



9A = 1-0 





Mod. B-H 


B-H 


ref [15 1 


ref |11] 


ref |28| 




0.084 


-0.051 






-0.062 




0.011 


-0.008 


-0.012 


-0.020/-0.035 





TABLE V. Results for Af£f 1 ' and M9£ m in MeV~ 
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FIGURES 
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FIG. 1. a) Diagrams included in the OEM approach, b) General diagram for the Standard and 
SU (4) calculation, c) Example of exchange diagram which is not included in the OEM. 
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FIG. 2. Calculated Gamow- Teller strength for the 48 Ca ^ 48 5c transition using the modi- 
fied Kuo-Brown potential (case a), the modified Bertsh-Hamamoto potential (case b), and the 
Bertsh-Hamamoto potential (case c). Each histogram represents the total strength in I MeV 
region versus the energy. 
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FIG. 3. Calculated (3 + strength for the 48 Ti -> 48 Sc transition using the modified Kuo-Brown 
potential (case a), the modified Bertsh-Hamamoto potential (case b), and the Bertsh-Hamamoto 
potential (case c). Each histogram represents the total strength in 1 MeV region versus the energy. 
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FIG. 4. Double beta amplitude versus the spin-orbit potential. All curves correspond to the 
modified B-H potential, except the Mg^- with squares which corresponds to the B-H potential. 
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FIG. 5. Detail of the contributions to the double beta amplitude, Mgq-, using the modi- 
fied Kuo-Brown potential (case a), the modified Bertsh-Hamamoto potential (case b), and the 
Bertsh-Hamamoto potential (case c). Each histogram represents the total contribution in 1 MeV 
region versus the energy. 
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